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ABSTRACT 

A comparison of polar solar wind proton flux upper limits derived 
using Saito's coronal density model, with Ly a measurements of the 
length of the neutral H tail of comet Bennet at high latitudes, shows 
that either extended heating beyond 2 R & is necessary some of the time 
or that Saito's polar densities are too low. Whichever possibility is 
the case, the fact that the solar wind particle flux does not appear to 
decrease with increasing latitude, indicates that the heavy element 
content of the high latitude wind may be similar to that observed in 
the ecliptic. It is then shown that solar wind heavy ion observations 
at high latitudes allow a determination of the electron temperature at 
heights which bracket the nominal location of the coronal temperature 
maximum thus providing information concerning the magnitude and extent of 
mechanical dissipation in the intermediate corona. 
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1. Introduction 

There is, at present, a lack of Information on the physical 
conditions in the polar regions of the solar corona and solar wind. This 
lack results in a corresponding uncertainty in the global characteristics 
and extent of that plasma of solar origin which fills interplanetary 
space and thereby controls the near solar environment. For example, only- 
little is known about the variation with heliographic latitude of so 
fundamental a flow parameter as the solar wind mass flux. Similarly, 
hardly any. ing is known about latitude -variations of the solar wind 
energy and m entum fluxes. Yet, these parameters may be very important 
in determining .he physical state of the polar corona and the size of the 
solar dominated cavity over the poles which separates the sun from the 
local interstellar medium. In addition, the existence of heavy elements 
in the polar solar wind may depend (Alloucherie, 1967, 1970; Geiss et 
al. , 1970) on whether or not the proton particle flux exceeds a 
temperature dependent lover limit. 

It is therefore useful to consider hypothetical variations of solar 
wind particle and energy fluxes with heliographic latitude. This task is 
approached in section 2 of this paper by calculating upper limit values 
of the polar solar wind particle flux implied by the most comprehensive 
coronal density model developed to date (Saito, 1970). This model was 
determined from an average K corona brightness distribution constructed 
using 15 solar eclipse observations as well as K-coronameter measurements 
all made at the minimum phases of the solar activity cycle. As a 
necessary result of the method employed, the model densities (and hence 
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the upper limit values of the solar wind particle flux derived below) 
determined for the polar regions are uncertain because it is not ^c=?ible 
to uniquely invert the convolution integral which relates the coronal 
brightness distribution to the average line of Bight electron density. 
Nevertheless, it is shown that if the polar coronal densities are as low 
as calculated using Salto’s model, then without extended heating, the 
emitted polar particle flux should be substantially less than that 
observed in the ecliptic plane at 1 AU and less than that necessary to 
drag coronal heavy elements away from the sun. However, limited evidence 
based on Ly a measurements of the neutral hydrogen tail of comet Bennet 
(Bertaux et al., 1973; Keller, 1973), is consistent with a polar solar 
wind flux at leaut as large as that observed in the ecliptic at 1 AU. 
These observations therefore require either an extended coronal heat 
source distinct from electron heat conduction or that Saito's polar 
densities are too low. In any event since the 

particle flux in the polar wind may be comparable to that observed in the 
equatorial wind, it is possible that coronal heavy ions at polar 
latitudes do indeed expand with the protons into interplanetary space. 

Since it is reasonable to expect that heavy elements will be 
observable in the polar solar wind, the range of ionization state 
"freezing in" distances is estimated In section 3 for selected heavy ion 
species at polar latitudes. It is found that the polar coronal density 
may be sufficiently low that the ionization states "freeze in" below the 
nominal location of the temperature maximum. Hence high latitude heavy 
ion observations may allow a determination of the thermal state of the 
intermediate and low corona and provide an estimate of the magnitude and 



extent of mechanical dissipation. Section h summarizes the main 
conclusions. 

2. Latitude Variations of the Solar Wind Particle Flux 

It is currently thought that the solar wind evolves from open field 
regions in the corui.a (see Hundhausen, 1972 for a review). Such regions 
are generally distinct from regions of activity and are generally characterized 
by low density. For th ». regions, the electron density, N, as a function 
of solar distance, r, and heliographic latitude, 0, has been modeled by Saito 
(1970) with the relation 



where R *» r/R and R is the solar radius. 

a » 

Upper limit values for the polar solar wind particle flux can be 
derived using relation 1 if various subsets of several reasonable 
assumptions concerning the state of the intermediate corona are adopted. 
These assumptions are l) the coronal gas consists of H, He and electrons 
only, 2) there is no extended heating other than that due to 

electron heat conduction much beyond the coronal temperature maximum, 

3) the energy equation may be closed with the standard relation 
Q * (Chapman, 195^; Spitzer, 1956) which assumes that binary 

coulomb interactions limit the mean scattering length of a thermal 
electron, h ) coronal electron and proton velocity distributions are very 
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nearly Maxwellian, 5/ wave-particle interactions and macroscopic wave 
pressure effects are negligible above the heating region, and 6) the 
magnetic field is open but not necesearily radial. 

The purpose of this section is to show that if -he coronal'density 
over the pole drops off as quickly as implied by Saito's analysis then 
some of the above assumptions may be tested by in situ solar wind 
observations. We begin with a standard single fluid formulation of t>s 
coronal expansion using equation 1 in place of an energy equation and 
derive upper limits for the particle flux at 1 AU. A separate treatment 
based on various possible forms of the energy equation is considered next 
to provide independent estimptes of the 1 AU flux upper limit. The 
results of this analysis are in agreement with those obtained by Durney 
and Hundhausen {197**). As will be shown in section 3 these upper limits 
are substantially lower than that observed in the ecliptic at 1 AU and 
are sufficiently low, that if all of the above assumptions are correct, 

He and many of the heavier ions should not expand with the protons away 
from the sun at polar latitudes. 


(i) Mass Flux. Momentum Flux and Density Equations 

The mass and momentum conservation equations are respectively; 

HVACR) = F 


"p""! = ( “ kT) 


NGM m M 
e p 



( 2 ) 

(3) 


2 

Here A(R) is the area of a flux tube which varies as R if the expansion is radial, 
G is the gravitational constant, is the mas3 of the sun, m^ is the proton 
mass, M is the mean molecular weight » (1 + 4a)/ (2 + l>a) where a is the He 
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abundance by number, R is Boltzmann's constant, H is the proton density, V is 
the bulk convection speed and T is the one fluid temperature. Concentrating on the 


region in the intermediate corona between R » r/R = 2 and 4, equation 

© 

1 for 6 = 90° can be simplified to the form: 


H(R) a 


7.9 x 10 -3 

2 cm 


If it is assumed that A(R) varies as IT then equations 2, 3, and U can be 


integrated analytically to obtain t(r) 



Here the subscript o refers to parameters evaluated at the base radius R q , 

In the following R is chosen equal to 2. 

© 

Eauation 5 can be rewritten in simplified form as follows: 

T(») = (|-) 6 (T - C [1- (Ja ) 7 ]- c 2 (nv) 0 2 t(|-) 6 - 2S - ill M> 

O O 

Here and are constants which are readily evaluated by comparing 
equations 6 and 5. Inspection of equation 6 shows that T(R) depends parametrically on 
two variables, T q and (NV) q . Following the analysis of Brandt et al. (1965) 
it is possible to show that two physically reasonable assumptions imply 
stringent constraints on the range of realizable values of T q and (NV) q . These 
two assumptions are: l) the derived temperature, T(R), must remain positive through- 

out the range of validity of equation 1; according to Saito (1970), R < U, 2) there 
is not sufficient external heating > /ond R = 2 to produce a second peak in T(R). 
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It is seen from the third tern on the right hand side of equation 5 
that for a constant T , increasing (NV) q eventually drives T(R) negative. 
The radius at which this happens can be increased beyond R = It by in- 
creasing T q . However if T q is too large, the (R/l^)^ term in front on the 
right hand side produces a second peak in T{R) beyond R = 2. Therefore 

acceptable ranges of T and (NV) can be determined as follows. The 

o o 

minimum value of T q , T^, is calculated for (NV) q = 0 under the assumption 
that T(r) < 0 for R > R^ where R^ is the limiting distance of validity of 
equation 1. This gives 



(7) 


Given R = ?, M = 0.5**7 (corresponding to a t!? lie abundance by number) and 
o 

assuming R„ = 3 and 3.5 then T_ = O .85 x 10^K and O .89 x 10^K respectively. 

A u 

Upper limits for T q arid r :.vi Q ere determined from equation 5 by finding 

the largest value of T q and (HV) q such mat T{R) > 0 and dT/dR < 0 for R in the 

range R < R < R„. Thus for each R the following two relations must be satisfied 
0 o X 



and 



Inspection of equations 8 and 9 shows that for (NV) q = 0 both conditions 

can be satisfied simultaneously. However for each R both conditions cannot 

be satisfied if (NV) is larger than some maximum value. This maximum is 

o 

obtained by equating the right hand sides of equations 8 and 9. 


N 


(NV). 


/GM \ 1/2 f /R \ 7 /R \ 7 
. „ % r/R..\6-2S t, \ , 


1/2 


{»)Rr-W(tTT 


2 <S <6 (10) 


p 

Setting S = 2 (radial flow) and (NV) = (NV) (R /R ) (the subscipt e refers to 

e o o e 

parameters evaluated at the orbit of the earth) , equation 10 is plotted 
in Figure 1 for R x * 3 and 3.5. The minimum value of the right 
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hand side of equation 10 for R in the range R q < R < jr^ is the maximum flux 
at 1 AU consistent with the assumptions dT/dR < 0 for R > 2 and T > 0 for 

Op -I f. 

R < R r Thus for R x = 3.0, (NV) g < 0.U x 10 cm sec with 0.8^ x 10 K < 

T(R = 2) < 1.1 x 10 6 K and for R =3.5, (lIV) < 0.2 x 10 8 cm^sec" 1 with 

e 

0.89 x 10 K ^ T(2) ^ O.98 x 10 K. The curves for T(R)^ corresponding to 
values of (NV) g and T(2) determined from equations 9 and 10 evaluated near 
the minimum of the curves in Figure i, are drawn in Figure 2, Drawn also 

for comparison are the polar scale height temperature , T^CR) = (GM m M)/(kr dinlT/dr) 

-2/7 P 

and the curve T a R ' . 

_2 

Since it is likely that the flow is more divergent than R inside of 

seme radius, Rq, it is necessary to consider how this possibility affects. 

the upper limit of (NV)^. This may he accomplished by assuming the area 

S 2 

of a flux tube increases as R out to and then as R from there to 1 AU, 

R . Using this model, (NV)_ = (NV) (R /? ) S = (NV) (R /R ) 2 and hence 
e D o o d e e p 

(NV) can be determined from equation 9 using the relation 


(NV) = 
e 



R <R <R < R 
0 X D 


( 11 ) 


Investigations of equations 10 and 11 for S in the range 2 < S < U, 

= 2, = 1.75 and R q < R < R x show that the maximum flux at 

1 AU consistent with a single temperature maximum below r is not significantly 
changed from its value for S = 2 ((NV) e < 0.2 x 10 cm sec x ). However, if 
S is sufficiently large and/or (R^/R ) is sufficiently small, this upper 
limit is raised. For example choosing S = 5 with ( R e/R 0 ) * 2 (which is 
equivalent to expansion from a polar region defined by 60° < 0 < 90° at 
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0 2 *■*! 

R q to the full hemisphere at y, CNV) e < 0.28 x 10 cm" sec" . It should 
be noted though that for all cases of nonradial expansion (NV) q is significantly- 
raised over that obtained for S = 2. 

ii) Energy Flux Supply to the Polar Wind 
An alternative approach to the polar particle flux problem is possible by- 
considering the energy equation. Here, a limit on (NV) g may be established if the 
velocity at 1 AU is known and if the usual assumptions about the state of the 
intermediate corona are made. Using a one-fluid, steady-state, spherically 
symmetric model, the energy equation 

iS)4B WB) «i-- m l + l' aV S U2 > 


may be combined with equations 2 and 3 and integrated to yield 


A(R) Q + F 


m MV 
P 


fkT 


GM m Ml 

®_£_ 

R R 

® -1 


e = constant . 
o 


Cl3) 


If a supersonic solar wind exists at 1 AU but not at F. then the dominant term at 1 

GM m M 

AU is F [4 m MV 2 ] while at R , the dominant terms are A(R )Q + F [t- kT ]. 

2 p e o o d c ri k 

r o o 

Therefore 


\ m MV 2 
2 p e 


Q o 

(NV) 



GM ra M 
9.R- 
R R 


o o 


(lM 


Further progress is not ’possible without an additional closure 
relation which gives Q q in terms of the lower velocity moments. Usually 
the Spitzer conductivity is assumed valid so that 

— 5/2 — 

Q = - K T 7/ VT 
H o 
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with k q = 7.7 x 10~ 7 erg cm' 1 sec' 1 K~ 7 ^ 2 (Chapman, 195^; Spitzer, 1956). 

However, it is also possible that the density is sufficiently low over 

the poles that equation 15 is not obeyed. In particular, it is possible 

that the polar density is so low that the dimensionless third moment, 

q = Q/[l.5 NkT(kT/m e ) ' ], becomes impossibly large at a low altitude 

(Parker, 1961*). For example if N q = 1.23 x lO'* cm ^ (Saito, 1970), 

T q = 0.98 x 10 (see section (i) above) and T K R 2 ^ 7 then q = 

(0.15)(R/R or q>l when R > 1.5 R • It is therefore probable 

o o 

that below this altitude instabilities develop (Forslund, 1970) which 

will limit Q to a value less than the Spitzer upper limit. In other 

words, the heat flux will be limited within 1.5 R Q thus effectively 

producing an isothermal region at lower altitudes and a region of steeper 
-2/7 

than R temperature decrease at higher altitudes. 

It is thus not clear how to estimate the value of Q in equation 

o 

lU. For the sake of concreteness two alternate approaches are adopted 

-2/7 

below. The first assumes equation 15 to be valid with T * R and the 
second adopts an exospheric approach. In both cases the solar wind Ke 
abundance is assumed to be a free parameter since its value is observed 
to be highly variable in the ecliptic at 1 AU and is not known at high 
polar latitudes. Such an assumption is necessary since, in contrast to 
the analysis presented in section (i) where upper limits for (NV) q were 
independent of M (see e.g. equation 10), the magnitude of the He 
abundance may be significant here. This fact results because most of the 
energy needed to drive the solar wind expansion goes into gravitational 
potential and kinetic energy which are both mass dependent. However this 
effect is more than compensated for by the fact that maximum values of 

-118- 



T q derived from the analysis in section (i) scale linearly with M (see e.g. 
equation 9). 

-2/7 

Assuming first that equation 1; is valid, T « R and 

T o =(0.98 x 10^)(M/0.5^7) K (see Figures 1 and 2) then upper limit values 

for (NV) can he calculated from equation lk for chosen values of M and 
e 

V £ . The results are summarized in Table 1 under the label (HV) e (Spitzer) 

for V g * 320, U50 and 750 km/sec and '! values corresponding to a He 

abundance of 0, 0.0U and 0.08 by number. 

In the appendix, an analysis is presented which shows that it is not 

clear whether or not electrons are ccllisionless below R = R^. If indeed 

coronal electrons are collisionless near to but outside of 2 R then Q 

e 

must be calculated u.3ing exospheric theory (Jockers, 1970; Lemaire and 
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Scherer, 1971a ,b; “Schulz and Eviatar, 1972; Hollweg, 197^; Eviatar and 
Schulz, 1975). In this approach, only those electrons above the electric 
potential harrier, |eA<J>| vith velocities directed away from the sun can 
carry heat. The three fluid energy equations may be combined and integrated 
to yield 


Pv)] 


(16) 


where the A symbol signifies a difference between any two radial distances 

and e is the electronic charge. Choosing R and R = 1 AU as the two 

o e 

ref-erenee-distancca then 


[ , „ 04 m 1 

o ® J 


(17) 


If both V g and the shape of the electron distribution at R , f(V), are known 
then Q q is readily evaluated using the relation 


oo Tf/2 TT/2 

%■/ / / (|-m e V 3 cos0 cos<J>)f(V)V 2 dV cos9d0d<|> (18) 

-n/2 -ir/2 


where — m e Vg 2 = |eA(J>j. Assuming a Maxwellian shape for f(V) then 



( 39 ) 
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Equations lU, 17, and 19 can be combined to give a self-consistency 
condition for 3 * Q q /[ (NV) Q kT o ] and hence an upper limit for (NV)" if 
the bulk convection speed at 1 AU is to be greater than or equal to V , 


Cw) o 



-CS+2.5) 

e 


[g + 7 + 


13.25 , 

e 1 


where 


, Mu V 2 GM 

- i P e - + 

2 kT RRkT 

o o © o 


Mm 

©j 


J and V > 0 
2 e 


( 20 ) 


( 21 ) 


Using equations 20 and 21 along with the assumptions that R q = 2, N = 1.23 

x 10^ cnf^ (Saito, 1970) and T q = 0.°8 x 10^(M/0.5^7) K (see e.g. Figures 1 

and 2), upper limit values for (NV) g have teen calculated for various values 

of M and V and are also listed in Table 1. 
e 

A comparison of the exospheric upper limits with the Spitzer upper 
limits for (NV) e shows that if both the wave-particle collision frequency and 
the coronal electron density are low enough over the solar pole so that an 
exospheric formalism is appropriate, very severe upper limits can be placed on 
the solar wind flux at 1 AU whether or not the He particles expand with the 
plasma. These upper limits fall well below that calculated from the mass and 
momentum equations alone. However, since it js likely that the corona is 
sufficiently turbulent that an exospheric formalism is not appropriate, the 
true upper limits for (NV) g may be less than but closer to that calculated 
using the Spitzer conductivity, 

A comparison of the upper limit values for (NV) e derived using the form 
of the energy equation which incorporates the Spitzer conductivity, with that 



derived using the mass and momentum equations* requires a knowledge of M and V . 
Reasonable choices for values of these quantities are made as follows. First, 
inspection of Table 1 shows that by choosing the solar wind He abundance to be 
h% by number, the error made in estimating (KV) e is probably less than 20 %. 
Therefore, for the purposes of this comparison, M is chosen to be 0.5**7. Con- 
cerning the speed of the polar solar wind at 1 AU several pieces of evidence 
have recently indicated that V g over the pole is higher than that observed in 
the ecliptic (Cole, 197M and may be close to 750 km/sec (Gosling et al., 1976; 

Feldman et al., 1976). If this is the case then from Tabl' ", energy considers- 

8 o 1 ^ 

tions require that (NV) be less than approximately 0.1* x 10 cm”"" sec - . This 

e 

value compares favorably with that derived using the ru.ss and momentum equations 

Op *1 

((NV) £ (0.2 to 0.1*) x 10 cm” sec” ). It is therefore concluded that if heat 

conduction is the dominant node of energy transport at about 2 R^, and if 

Saito's polar density model is correct then the particle flux of the polar 

8 2 1 

solar wind should be less than about 0.5 x 10 cm” sec” . This upper limit is 
about a factor of 7 times less than the solar wind particle flux observed in the 
ecliptic at 1 AU (Feldman et al., 1976). 

3. Latitude Variations of Heavy Ion "Freezing In" Distances 

Alloucherie (1967) derived an approximate criterion necessary for a 
heavy ion of mass Am^ and charge Z to diffuse upward in an expanding corona. 


It should be noted that this value is an upper limit. If the heat flux, Q, 

is regulated below 1,5 R as suggested earlier, then the region below 1.5 R 

o o 

becomes more nearly isothermal thereby reducing VT, Q, and hence the upper 

limit value derived for ( NV ) . 

e 


- 122 - 



His result is: 


(NV) q > 



B/n e lnX(Z^/A) 


( 22 ) 


where the symbols are as previously defined and lnX is the coulomb 

r n 6 

logarithm. Choosing N q = 1.23 x 10 5 cm , T q * O .98 x 10 K and expressing 
22 in terms of (NV) e we get: 


(nv) > 

e 


o 

2*0 x 10 
(Z 2 /A) 



-2 -1 


cm sec 


(23) 


Since for a radial expansion (S = 2) this limit is approximately a 
factor of U to 10 times greater than the upper limit for (NV)^ derived above, 
it is reasonable to conclude that if Saito's model is correct and if the polar 
corona is not externally heated above r = 2R ffi , He may not expand with the 
solar wind. This conclusion remains ^alid for S < 1 and (R^/R^) = 2 as veil. 

However, observation of the length of the neutral hydrogen tail of Comet 

Bennet (Bertaux et al., 1973; Keller, 1973) as a function of heliographic 

latitude indicates that the polar solar wind flux is at least as large as 
8 —2 — 1 

2 x 10 cm sec at 1 AU. This value is in disagreement with the upper 
limits deduced in section 2. It is therefore concluded that at least one 
of the assumption., made in the above analysis is not correct and that it 
should indeed be possible to observe solar wind heavy ions at polar latitudes 
at 1 AU. If true then measurements of the populati m densities of individual 
heavy ion ionization states will yield information concerning the temperature 


- 123 - 



structure of that region in the polar solar corona where the various ionization 
states "freeze in." 

It is possible to determine the "freezing in" distances of the various 

heavy ion species as a function of heliographic latitude using equation 1 if 

the following assumptions are made: l) the flow is radial; 2) the velocity 

distribution is Maxwellian; and 3) the electron temperature, T, depends on the 

radius, r, as T * T (r/R )“ Y . Following previous work 
$ © 

(Hundhausen et al., 1968a, b; Bam> jt al., 1974) these distances are 

defined as those for which the expansion rate, T * (YdlnN/dr) becomes 

larger than the ionization state changing rate, + T c J «• N(R^ + C^). 

He-: ■’ is the electron temperature at the base of the corona, V is the 

solar wind speed, is the rate of recombination from state 1 to state i-1 

and is the rate of collisions! ionization from state i to state i + 1.* 

Changes in the "freezing in" distances with latitude of a sample of the 

most abundant ions are shown schematically in Figure 3 superimposed on scale 

height temperatures calculated using equation 1 for 9*0° and 90°. For 

purposes of illustration an isothermal corona with T * 1.0 x 10^K and a 1 

r 

AU particle flux of 2.5 x 10 cm sec were assumed for evaluating + x~ . 

The scale height temperature for a static corona is given by T * 

a 

2 

(GM m M)/(kr dlnN/dr). Inspection of Figure 3 shows that the region in the 
o p 

corona for which temperature values can be determined frc.a solar wind heavy 
Ion data moves inward from above to below the temperature maximum as 9 
varies between 0° and 90°. Thus at some intermediate latitude, coronal 
temperatures bracketing the maximum can be sampled allowing the magnitude 
and extent of mechanical dissipation in the Intermediate corona to be 
estimated (see e.g. the analysis of Brandt et al., 1965). 

*Colli8lonal ionization and radiative recombination (including dielectronic 
recombination) coefficients for 0, Si, and Fe were kindly supplied by Dr. A, Dupree. 



1». Summary and Conclusions 

In this paper two related aspects of the physical state of the inter- 
planetary plasma at high solar latitudes were explored. In the first part 
upper limits for the polar solar wind particle flux were derived using a set 
of reasonable assumptions concerning the base coronal conditions along with 
Saito's (1970) coronal density model. In the second part, it was determined 
whether this flux was sufficient to drag the heavier ions away from the sun 
into interplanetary space. 

From the analysis in the first part it was concluded that if Saito's 
model is correct, the polar electron density is sufficiently low that in the 
absence of extended heating the solar wind flux at high latitudes should be 
at least a factor of from h to 10 tines less than that observed in the 
ecliptic at 1 AU. Such a low particle flux was shown in the second part 
to be small enough that most heavy ions would not be expected to expand 
with the protons into interplanetary space. 

However, indirect and limited evidence available at present is con- 
sistent with a polar solar wind that Has at least as large a velocity (Coles 
et al., 197k; Brandt et al., 197M and as large a particle flux (Bertaux 
et al., 1973", Keller, 1973) as that observed in the ecliptic at 1 AU. From 
the analysis presented in section 2, these observations then require either 
that extended heating distinct from that provided by electron heat conduction 
is necessary some of the time above 2.3^ or that Saito's polar densities 
are too low. Whichever is the case, the fact that the 

solar wind particle flux does not appear to decrease with increasing heliographic 
latitude (Bertaux et al., 1973; Keller, 1973) indicates that coronal heavy 
ions may be expected to expand with the protons away from the sun. If true 
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then measurements of the population densities of individual heavy element 
ionization states in the polar wind vill provide information at 1 AU con- 
cerning the thermal state of that region in the intermediate corona where 
the respective ionization states freeze in. It turns out that the latitude 
variation of these freezing in distances calculated using Saito's model 
is such that the region in the corona for which temperature values can be 
determined moves inward from above to below the nominal location of the 
temperature maximum as 0 varies between 0° and 90°. Therefore, measurements 
of heavy ions at high solar latitudes may provide valuable information con- 
cerning the magnitude and extent of mechanical dissipation in the intermediate 
polar corona. 
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Table 1 

Upper Lir.it Values of (HV) e 
Consistent With the Enerry Equation 


km/sec 

% 

-2 -1 

cm sec 

-2 -1 

cm sec 

V 

e 

He/K 

(NV) (Spitzer) 
e 

(NV) Exospheric 
e 

320 

0 

1.02 x 10® 

1.2U x 10 6 

U50 

0 

O 

0.70 x 10 

U.96 x 10** 

750 

0 

0.33 x 10 8 

0.U1 

320 

k 

1.28 x 10 8 

c 

1.3 x 10 

U50 

u 

0.88 x 10 8 

u 

5.2 x 10 

750 

h 

0.U1 X 10 8 

o.u 

320 

8 

O 

1.5 ! ‘ x 10 

l.fco x 10 6 

U50 

8 

O 

1.06 x 10 

5.59 x 10 U 

750 

8 

O 

0.1*9 x 10 

0.1*6 
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Appendix 

Comparison of Expected Electron-Electron Collision Lengths 
with Scale Lengths in the Polar Corona 
The magnitude of the electron conductivity in the polar corona 
depends critically on the electron-electron collision length,. If S- c 

is small enough then the Spitzer conductivity is applicable but if it is 
too large, then an exospheric approach is needed to evaluate the polar 
electron heat flux. It turns out that, according to Saito, H q is 
sufficiently low over the pole that it is not clear vhether or not * v .ermal 

electrons are collisionless above R = 2R , For example the self scattering time for a 

o 

thermal electron at 2R^ is T c = (l.l x 10 -2 )T /N = 87 sec (Spitzer, 1956) 
whereas at that distance the expansion time (assuming a radial magnetic 
field) is T g = [ (kT o /m e )^ 2 dlnN/dR] ^ = 60.5 sec. Furthermore the coulomb 
scattering length, J , defined ^y 


R + l 
o e 


1 = 


(kT /m ) 1/2 
o e 


/ 


dR 

T 


(Al) 


may be either larger than or smaller than the temperature scale length, 

= I-dlnT/dR) = 3-5 R Q , depending on the value of the maximum altitude 
at which Saito's density model is valid, R x> This is readily shown by 
assuming a density model consistent with Saito's results (1970): 


k-6 


» ■ \ (fj 

-.©'‘tr 


2 < R < R„ 


R x < R 


(A 2) 
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with Combining A1 and A2 and using R q = 2 with N q from equation 1 

evaluated at 8 = 90° it is found that A c /R 0 ■ 1-33, 2.U8 and 6.69 for R^ = 
3-0, 3.5 and U.O respectively. 

Since the actual value of Z may be either less than or greater 
than depending on the value of R^, it is not known whether the Spitzer 
conductivity or a conductivity calculated using exospheric theory is most 


valid in the polar solar corona. However, the fact that & c is of the same 
order of magnitude as suggests that neither of the above is correct and 
that to obtain an accurate determination of the true conductivity a 


kinetic approach may be necessary. 
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Figure 1. 


Figure 2 . 


Figure 3. 


Figure Captions 


Plots of equation 9 for R^ = 3 and 3.5* Radial flow is assume^ 

p 

and hence (NV) * (NV) (R /R ) . The minimum value of each 
e o o e 

curve corresponds to the maximum flux at 1 AU consistent with 
the assumptions dT/dR < 0 for R 2 and T > 0 for R ^ R^, 

Coronal temperatures, T(R) for values of (NV) g and T(2) determined 
from equations 9 and 10 evaluated at the minima of the curves in 
figure 2. Drawn also for comparison are the polar scale height 
temperature, T„(R) and the curve T a R ' . 

n. 


Variations with, latitude of the ''freezing in" distances of a 
sample of the most abundant heavy ions. Scale height tempera- 
tures are calculated using equation 1 for 0=0° and 90 ° end a 
constant "freezing in" temperature of T^ = 1 x 10^K is assumed. 
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